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Abstract. We report high-sensitivity millimetre observations of several molecular species (^^CO, HCN, HNC, CN, 
HCO^ and N2ll^) in a sample of compact planetary nebulae. Some species such as HCO^ and CN are particularly 
abundant compared to envelopes around AGB stars or even interstellar clouds. We have estimated the following 
average values for the column densities ratios: CN/HCN~2.6, HCO+/HCN~0.5, and HNC/HCN~0.4. Thus, the 
chemical composition of the molecular envelopes in these compact PNe appears somewhat intermediate between 
the composition of proto-PNe (such as CRL2688 or CRL618) and well evolved PNe (such as the Ring, M4-9, 
or the Helix). From observations of the CO isotopomers, we have estimated that the ^^C/^"^C ratio is in the 
range 10 ^ ^^C/"C ^ 40. These values are below those expected from standard asymptotic giant branch models 
and suggest non-standard mixing processes. The observed molecular abundances are compared to very recent 
modelling work, and we conclude that the observations are well explained, in general terms, by time-dependent 
gas-phase chemical models in which the ionization rate is enhanced by several orders of magnitude with respect 
to the average interstellar value. Thus, our observations confirm that the chemistry in the neutral shells of PNe 
is essentially governed by the high energy radiation from the hot central stars. The complexity of the chemical 
processes is increased by numerous factors linked to the properties of the central star and the geometry and degree 
of dumpiness of the envelope. Several aspects of the PN chemistry that remains to be understood are discussed 
within the frame of the available chemical models. 
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1. Introduction 

While planetary nebulae (PNe) are primarily identified 
from their optical emission arising in the ionized gas, they 
I often contain a substantial amount of neutral gas which 
is a more direct relic of the mass lost on the asymptotic 
giant branch (AGB). CO and H2 surveys of PNe have 
revealed that a large fraction of the gas in massive PNe 
is in molecular form (Huggins et al. 1996, Kastner et al. 
1996). This molecular component may be dominant during 
a large part of the evolution of the nebulae, and thus plays 
a prominent role in the mass distribution and the shaping 
of PNe. 

The molecular gas in PNe is affected by the large ra- 
diation fields from the central stars which are expected 
to cause a dramatic effect on the chemical composition. 
The studies of the chemistry of the neutral gas in PNe 
have been very scarce so far. Cox et al. (1992) observed 
the evolved PNe NGC 6072 and IC 4406 at milHmetre 
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wavelengths, and detected HCN, HNC, CN and HCO+ in 
both objects. More recently, Bachiller et al. (1997a) made 
a more complete study, covering an evolutionary sequence. 
If we except the two PPNe of their sample, five PNe were 
searched and detected over this large variety of molecular 
species: the young PN NGC 7027 and the evolved ones 
NGC 6720 (the Ring), M 4-9, NGC 6781 and NGC 7293 
(the Helix). The molecular content of proto-PNe has also 
been investigated by observing some selected objects such 
as CRL 618 (Bujarrabal et al. 1988 ; Cernicharo et al. 
2001), and OH 231.8-1-4.8 (Sanchez Contreras et al. 1997). 
It clearly appears from these studies that the chemical 
evolution is dominated by photodissociation, ion/radical- 
molecular reactions, and shocks. In particular, the sur- 
vival of some molecules at well advanced stages of the 
evolution is made possible thanks to the shielding of the 
radiation, which is directly connected with the presence 
of inhomogeneities in the envelope. Very dense clumps, 
directly observed as cometary globules in the nearest PNe 
such as the Helix (Huggins et al. 1992, 2002), provide an 
environment which is well suited to the protection of the 
molecules and of its on-going chemistry. The chemistry 
of the PNe neutral envelopes has been modelled by sev- 
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eral authors. The first work was done by Black (1978) 
who studied the neutral-ionized transition zone in an ho- 
mogeneous nebula. Howe et al. (1992, 1994) modelled 
the evolution of a dense clump similar to those of the 
Helix nebula. More recently, Hasegawa et al. (2000) and 
Hasegawa & Kwok (2001) reported detailed models of the 
neutral envelope of NGC 7027, and Ah et al. (2001) con- 
structed time-dependent gas-phase models which repro- 
duce well the molecular abundances in evolved PNe such 
as NGC 6781, M4-9 and the Helix, observed by Bachiller 
et al. (1997a). 

With the remarkable exception of NGC 7027, the pre- 
vious chemical studies carried out so far have been mainly 
devoted either to rather evolved PNe (such as the Ring 
or the Helix) or to proto-PNe (such as CRL618 or 
OH 231.8-1-4.8). We report here on high-sensitivity obser- 
vations of molecular lines in a sample of 7 compact, young 
or intermediate-aged PNe, which improve our knowledge 
of the chemistry during the nebular development. The 
sample and the observations are presented in Sect. 2. An 
account of the results is given in Sect. 3. These results are 
discussed in Sect. 4 making particular emphasis on the 
chemistry of N-bearing compounds, on the formation of 
molecular ions, and on the importance of the dumpiness. 



2. Sample and observations 

2.1. PNe sample 

Our sample is made of 6 compact PNe covering a wide 
range of evolutionary status and whose molecular chem- 
istry has not been studied up to now. These are BV 5-1 
(see Josselin & Bachiller 2001 for its CO content), K 3-94, 
K 3-24 (Josselin et al. 2000), M 1-13, Ml-17 and IC 5117 
(Huggins et al. 1996 and references therein). Although it 
was already observed by Bachiller et al. (1997a), as repre- 
sentative of young PNe, we also observed NGC 7027 for 
comparison and calibration purposes. Table ^ lists the 
main properties of the observed nebulae: the most usual 
name (Col. 1), standard galactic name (2), kinematic dis- 
tance - based on the simple approximation for the velocity 
structure of the Galactic disk of Burton (1974) - (3), op- 
tical size of the ionized nebula (4), radial and expansion 
velocities measured in the C0(2-l) line (5, 6), integrated 
intensity of the C0(2-l) line (7), mass of the molecu- 
lar envelope (8), and ratio Mm/Mj between the molecular 
mass and the ionized mass (9), which is a good indica- 
tor of the evolutionary stage, almost independent of the 
distance (Huggins et al. 1996). Concerning the distance, 
we have preferred the kinematic estimate over the statis- 
tical estimate, since we believe that the former is more 
precise in these cases (see Josselin et al. 2000 for details) . 
Nevertheless, for NGC 7027, the best studied object in 
the sample, in accordance with previous works, we have 
assumed the standard estimate of 0.7 kpc (see e.g. Huggins 
et al. 1996).In the case of M 1-13, M 1-17 and IC 5117 the 
values of Mm and Mm /Mi have been re-scaled from those 



given by Huggins et al. (1996) according to our choice of 
distance (Mm oc distance^ and Mm/Mi oc distance~^/^)Q 

2.2. Observations 

The observations were carried out with the IRAM 30-m 
radiotelescope on Pico Veleta (near Granada, Spain) in 
July 1999. The search for molecular emission was made 
simultaneously in four lines. We used SIS receivers oper- 
ating in the bands around A 3 mm and A 1 mm with typical 
system temperatures of 150 and 300 K, respectively. The 
half-power beam width of the telescope is 22" and 10" at 
3 and 1 mm respectively. The spectrometers were 256x1 
MHz filter banks providing a velocity resolution of 2.6 and 
1.3 km/s at 3 and 1mm, respectively. The observations 
were done by wobbling the secondary mirror to a distance 
of 2' from the source. The pointing accuracy, checked ev- 
ery hour, was within 3". Calibration was achieved with a 
chopper wheel. The intensities given hereafter are in units 
of main beam brightness temperature. 



3. Results 

3.1. Spectra 

Spectra obtained toward the sample of PNe are shown 
in Fig. 0. Table ^ lists the line integrated intensities and 
upper limits derived from such spectra. The integrated 
intensities for NGC 7027 are generally found to be com- 
patible with those previously reported by Bachiller et al. 
(1997a). A significant difference was only found for the 
^^C0(2-l) line, for which we report an intensity higher 
by a factor of 30%. We caution that slight pointing errors 
can result in important differences in the line intensities 
of such compact objects, and we believe that the observa- 
tions reported here were made with better pointing accu- 
racy. 

HCO+ is detected in all the PNe of our sample. At least 
one of the N-bearing molecules we looked for (HCN, HNC, 
CN) is also detected in every PN, the three molecules be- 
ing detected in K 3-94 and M 1-17. The HNC(l-O) emis- 
sion is weak, and we only have at best tentative detections. 
"CO is detected in BV 5-1, M 1-17, K 3-24, IC 5117 and 
NGC 7027. 

We also observed the N2H+(l-0) line in the sample. 
The first detection of this line in a circumstellar enve- 
lope was reported by Cox et al. (1993) in NGC 7027 (see 
their Table 1). A high-sensitivity spectrum taken toward 
NGC 7027 is shown in Fig. |[ and its intensity is found 
to be consistent with that reported by Cox et al. (1993). 
This line was not detected in any other PN. The detection 
limits are reported in Table 



the ionized mass is derived from the A6 cm emission and 
the angular nebular size, resulting in a distance^^^ dependance, 
as shown by Gathier (1987). 
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FIN (j 


Uk 
(kpc) 
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(km s ) 


\r 

* cxp 

(km s ) 


T 

-'-12C0{2-1) 

(K km s ) 


Mm 
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(10 Mq) 


Mm/Mi 


BV 5-1 


119.3-1-00.3 


5.5 


20 


-73 


10 


9.7 


2.8 


0.3 


K 3-94 


142.1-^03.4 


5.5 


7 


-69 


16 


14.1 


3.4 


0.2 


M 1-13 


232.4-01.8 


2.0 


10 


27 


18 


22.2 


25.4 


1.8 


M 1-17 


228.8-f05.3 


2.5 


3 


28 


39 


66.2 


6.1 


14.0 


K 3-24 


048.7-f02.3 


3.5 


6 


44 


24 


23.2 


2.1 


0.2 


IC 5117 


089.8-05.1 


3.0 


1 


-11 


14 


19.6 


24.5 


2.1 


NGC 7027 


084.9-03.4 


0.7 


14 


26 


23 


278.0 


22.0 


9.4 




Fig. 1. Spectra of the hnes observed in the PN sample. The temperature scale is the same for every line for each PN, 
except for the HCO+ line in PNe BV 5-1 and K 3-94, where the data have been expanded by a factor of 10. The 
velocity scale of the CN profiles are calculated with respect to the frequency of the strongest component of both lines 
(z^o = 113.490945 GHz for the CN(l-G) hne and = 226.874766 GHz for the CN(2-1) hne). Some of the data have 
been smoothed for display purposes in the present figure (line intensities reported in Table 2 have been measured at 
the original resolution). 



3.2. Column densities and their evolution 

Estimates of the column densities have been made under 
the usual assumptions of LTE and optically thin emission. 
The excitation temperature was assumed to be uniform 
and equal to 25 K, which is the kinetic temperature de- 
duced from CO observations in a few other nebulae (see 



Bachiller et al. 1993). When the line was not detected, we 
used three times the upper limits to the integrated intensi- 
ties given at a 1(7 level in Table |^ to derive upper limits in 
the corresponding column densities. The CN column den- 
sities have been calculated from the CN(2-1) line intensi- 
ties, since this line is expected to be optically thinner than 
the CN(l-O) line and, in addition, the 2-1 observations 
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Table 2. Integrated intensities of the observed lines, in main-beam scale. Intensities are in K km s ^ and 1 rms 
errors are given in parenthesis. The quoted values for the CN lines refer to the high frequency fine-structure groups 
J = 3/2 ^ 1/2 for the CN(l-O) line and J = 5/2 ^ 3/2 for the CN(2-1) line. 



Name 


^^CO(2-l) 


HCN(l-O) 


HNC 


1-0) 


CN(N= 


= 1-0) 


CN(N= 


=2-1) 


HCO + 


(1-0) 


N2H+(l-0) 


BV 5-1 


1.10 


(0.05) 


0.54 


(0.02) 




(0.03) 


0.77 


(0.03) 


2.22 


(0.08) 


0.12 


(0.02) 


(0.03) 


K 3-94 




(0.15) 


0.41 


(0.03) 


~0.35 


(0.09) 


0.96 


(0.05) 


3.54 


(0.10) 


0.18 


(0.02) 


(0.04) 


M 1-13 




(0.2) 


1.16 


(0.07) 


~0.30 


(0.07) 




(0.2) 


2.16 


(0.3) 


0.95 


(0.07) 


(0.05) 


M 1-17 


2.49 


(0.3) 


2.77 


(0.11) 


~0.36 


(0.11) 


1.39 


(0.2) 


2.77 


(0.6) 


1.18 


(0.08) 


(0.2) 


K 3-24 


~0.66 


(0.11) 


0.58 


(0.02) 




(0.02) 


0.61 


(0.05) 


3.58 


(0.12) 


0.42 


(0.02) 


(0.01) 


IC 5117 


0.85 


(0.11) 




(0.03) 


~0.15 


(0.04) 




(0.03) 


0.81 


(0.08) 


0.45 


(0.07) 


(0.06) 


NGC 7027 


26.30 


(0.2) 


12.56 


(0.2) 




(0.2) 


19.15 


(0.3) 


45.49 


(0.8) 


29.78 


(0.4) 


1.74 (0.06) 




-50 50 100 

velaclly (km/s) 



Fig. 2. Spectrum of the N2H+(l-0) line towards NGC 
7027. 



provide higher resolution. It should be pointed out that 
the hyperfine structure of the CN spectrum may introduce 
some additional uncertainties in the estimates of column 
densities. This is because of the possible anomalous exci- 
tation through near-infrared pumping which can produce 
weak masers (Bachiller et al. 1997b). Unfortunately, the 
separation of the different hyperfine lines is very difficult 
in the case of our nebulae because of the complex kine- 
matic structure of the line profiles, which often include 2 
or 3 components at close velocities (see for instance the 
CO profiles in Huggins et al. 1996 and Josselin et al 2000). 
We thus simply assumed that the hyperfine line compo- 
nents within a fine group presented standard LTE ratios. 
The resulting column density estimates for all observed 
molecules are given in Table | 

For all the observed species, in all objects other than 
NGC 7027, the column densities are relatively small which 
is an indication that the hypothesis of optically thin emis- 
sion is probably appropriate. Even in NGC 7027, the large 
column densities may not be incompatible with opti- 
cally thin emission (Hasegawa & Kwok 2001). However, 
the assumptions of LTE and uniform excitation tem- 
peratures are much more questionable. In fact, all the 
observed species but ^'^CO have high dipole moments 



(^ 1.5 Debye) and are thus expected to present sub- 
thermal populations. 

In order to study the evolution of the abundances in 
the different nebulae, we thus decided to use ratios with 
respect to HCN (which was detected in all nebulae but 
IC5117). Since HNC, CN and HCO+ have high dipole 
moments similar to that of HCN, non-LTE effects on the 
column densities should be of the same order of magnitude 
for these species so the column density ratios should be less 
affected by these uncertainties. Furthermore, Hasegawa 
& Kwok (2001) found that HCN, CN and HCO+ emis- 
sion regions have similar sizes in the case of NGC 7027, 
while CO emission is more spatially extended. Beam dilu- 
tion effects should thus be less important when compar- 
ing emission from the observed species but CO. Figure |^ 
displays the values of the HNC, CN and HCO+ column 
densities with respect to that of HCN. The column density 
ratios are plotted as a function of the molecular to ionized 
mass ratio (Mm/Mi). Data from Bachiller et al. (1997a) 
for CRL 2688, a typical post-AGB star and the youngest 
object in their sample, and NGC 7293 (Helix, their oldest 
PN) have been added, in order to stress the evolution of 
the abundances. Concerning CRL 2688, one should keep 
in mind that this object is too young for ionized gas to be 
present, so its abscissa is arbitrary. Furthermore, its HCN 
line is optically thick so the relative abundances should be 
considered as upper limits. 

Since, as mentioned before, Mm/Mj is a good indica- 
tor of the evolutionary state of the nebula, the abscissas 
in Fig. ^ can be considered as a time axis, and the lines 
connecting the values of the column density ratios should 
reflect the evolution of the abundance ratios. It is interest- 
ing to note some general trends in the abundance ratios. 
The most striking point is the dramatic increase in HCO"'' 
abundance relative to HCN from post-AGB stars to PNe, 
by more than two orders of magnitude. A second-order 
noticeable effect is the slight decrease of the HCO'^/HCN 
ratio from rather young PNe (Mm/Mi ^ 10) to older ones. 
The CN abundance relative to HCN is observed to increase 
by a factor of about 20 from the youngest objects in the 
sample (NGC 7027 and M 1-17) to the oldest ones (Mm/Mi 
^ 0.2). We remind here that the rather high value found 
for CRL 2688 must be attributed to an under-estimation 
of the HCN column density. The increase in the CN abun- 
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Table 3. Column densities (in cm^^) derived from observed line intensities given in Table ^ (see Sect. 4.1 for details). 



Name N(^=^CO) N(HCN) N(HNC) N(CN) N(HCO+) 



BV 5-1 


3.9 10" 


1.3 10^^ 


<2.4 10" 


3.0 10^^ 


1.7 10" 


K 3-94 


< 1.6 10" 


1.0 10^^ 


~ 9.1 10" 


4.7 10^2 


2.6 10" 


M 1-13 


<2.1 10" 


2.9 10^^ 


~ 7.9 10" 


2.9 10^^ 


1.4 10^^ 


M 1-17 


8.9 10" 


6.9 10^^ 


~ 9.4 10" 


3.7 10^^ 


1.7 10^^ 


K 3-24 


~2.4 10^^ 


1.4 10^^ 


<4.2 lO" 


4.8 10^^ 


6.0 10" 


IC 5117 


3.0 10" 


<2.2 10" 


~ 3.8 10" 


1.1 10^^ 


5.5 10" 


NGC 7027 


9.4 lO^'^ 


3.1 10" 


<9.3 10^^ 


6.1 10" 


4.3 10" 



dance occurs essentially during the early stages of the PN 
evolution, the CN/HCN ratio becomes rather constant af- 
ter then. The HNC/HCN column density ratio does not 
present a very clear tendency with age. But it should be 
pointed out that the HNC line was only tentatively de- 
tected in some of the objects, and that more sensitive HNC 
observations are needed before drawing any firm conclu- 
sion. Nevertheless, this ratio seems again to be relatively 
lower in the two youngest objects (HNC/HCN < 0.2) 
than in some more evolved ones (HNC/HCN up to ~ 1). 
These values are in the range of those previously observed 
in other nebulae (Bachiller et al. 1997a) and estimated 
theoretically (Ah et al. 2001). 

The examination of Fig. ^ also gives information about 
the evolutionary stages of the nebulae. For instance, M 1- 
17, which was previously suspected to be a young PN 
because of the presence of wide wings in its spectrum 
(Bachiller et al. 1991), presents a high Mm/Mi ratio. This, 
together with its molecular abundances (see Fig. ||), con- 
firms now that M 1-17 is a very young object, similar in 
many respects to NGC 7027. On the other hand, BV 5-1, 
K 3-24 and K 3-94 seem to be at a well advanced evolu- 
tionary stage, close to that of NGC 7293. 

A remarkable exception among the observed be- 
haviours is the case of IC 5117, the only nebula not de- 
tected the in HCN 1-0 line, and with rather weak CN 
emission (only the CN 2-1 line was detected), which seems 
typical of an oxygen-rich chemistry, whereas ionized gas 
in IC 5117 appears carbon-rich (C/0^1.3, Hyung et al. 
2001). A C/0 ratio >1 is consistent with the idea that 
massive PNe as those in our sample result from the evolu- 
tion of stars with intermediate masses (M <^ 3 Mq) and 
thus have experienced the 3^^*^ dredge-up, which makes 
them carbon-rich. However the determination of the car- 
bon abundance in IC 5117 is very uncertain. Furthermore, 
the existence of massive oxygen-rich PNe is not excluded, 
because of hot-bottom burning which could convert ^^C 
into ^^C and ^''N, decreasing the relative abundance of 
carbon (Renzini & Voli 1981). One could also think about 
a chemical stratification in the nebula, as the molecular 
gas corresponds to mass ejection prior to that traced by 
ionized gas. Additional observations are needed to under- 
stand this peculiar PN. 



3.3. Carbon isotopic ratio 

In the case of PNe, the observed I(i2co(2-l))/I("CO(2- 
1)) intensity ratio is expected to be a good indicator of 
the isotopic ^^C/^^C abundance ratio (see Palla et al. 
2000). In fact, the intensity ratio equals the abundance ra- 
tio when the following assumptions are accomplished : (i) 
the 12 CO (2-1) and "CO (2-1) rotational levels are ther- 
mahzed, (ii) the ^^00(2-1) and ^^€0(2-1) are not very 
optically thick, and (iii) the beam filling factor is similar 
for both lines. Since the dipole moment of CO is small 
0.11 Debye), the first assumption is very plausible. In ad- 
dition, saturation effects are expected to be negligible for 
most evolved PNe, since the I(CO(2-1))/I(CO(1-0)) ra- 
tio suggests optically thin emission (Huggins et al. 1996, 
Josselin et al. 2000). But this may not the case for the 
youngest PNe (such as NGC 7027 and M 1-17), for which 
the i^CO lines are probably optically thick; isotopic ratios 
for these objects should thus been considered as lower lim- 
its. 

The validity of the hypothesis of equal beam filling 
depends directly on the spatial distribution of the CO iso- 
topomers, and thus on selective processes of formation and 
destruction of the molecules. Two types of reactions may 
be important to determine the CO isotopomer ratio: selec- 
tive photodissociation, ^'^CO being more easily photodis- 
sociated, and chemical fractionation through the reaction 

12^Q _^13 fj+ _^ 13(jQ _^12 ^+ 

But these effects are generally thought to compensate each 
other (see Mamon et al. 1988 for a more detailed discus- 
sion) . In summary, we believe that for most of the PNe dis- 
cussed here, i^c/^C = I(i2CO(2-l))/I("CO(2-l)). The 
isotopomer ratios resulting from the observations are sum- 
marized in Table ^ 

We find ratio values in the range 
10<i2CO/i3CO<40. The ^^C/^^C ratio was re- 
cently measured in some additional PNe by Palla et al. 
(2000) and Balser et al. (2002), and the values reported 
by these authors are in general agreement with the ones 
reported here. These values are below those expected 
from standard asymptotic giant branch models and 
suggest non-standard mixing processes. 
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Table 4. Carbon isotopic ratios in the PNe of our sample. 



Name 




BV 5-1 


9 


K 3-94 


k 31 


M 1-13 


Z 37 


M 1-17 


27 


K 3-24 


35 


IC 5117 


23 


NGC 7027 


11 



with the models they quote. We thus confirm that non- 
standard mixing processes should thus be important only 
in low-mass stars, as suggested by Charbonnel (1995). 



< ■ 



10 



10 1 

Mm/Ml 



10 



4. Discussion 

Wc next discuss the chemical processes which must be 
relevant in the case of the molecules observed here by 
strengthening the differences between the conditions pre- 
vailing in PNc from those of the interstellar medium. 
Special attention is devoted to the dumpiness of the neu- 
tral FN envelopes since clumping is expected to play a 
key role in the preservation and evolution of the molecu- 
lar species. 



Fig. 3. Abundances of HNC, CN and HCO+ relative to 
HCN as a function of the molecular to ionized mass ra- 
tio. The ratio is plotted with a filled circle when both 
molecules are detected and with an empty triangle if one 
of the molecules is not detected, the orientation of the tri- 
angle indicating if the limit is a lower or an upper one. 
Data for CRL 2688 and NGC 7293, indicated by squares, 
are from Bachiller et al. (1997a). Horizontal dashed lines 
indicate the average of the relative abundances within our 
sample. 



We have estimated the mass of the progenitors of BV 
5-1 and K 3-94, by assuming relevant stellar parameters 
taken from the Acker et al. (1992) catalogue, and by fol- 
lowing the procedure described in Palla et al. (2000). For 
M 1-17, IC 5117 and NGC 7027, estimates are provided by 
Palla et al. (2000). For M 1-13 and K 3-24, the lack of the 
central star photometry prevents any realistic estimate to 
be done. Of the studied sample, BV 5-1 and K 3-94 ap- 
pear to have the most massive progenitors, with ~5.5 and 
^4.0 Mq, respectively (we note however that the estimate 
of the progenitor mass requires the previous estimate of 
the luminosity of the central star which in turn depends 
on the choice of the distance, which as discussed above re- 
mains quite uncertain). Then, when reported in Fig. 3 of 
Palla et al., the ^^C/^^C isotopic ratio appears consistent 



4.1. The survival of molecules in PNe 

The compact PNe discussed here sample the evolution 
of the chemistry from the younger objects to rather de- 
veloped PNe. An interesting result from our observations 
is the continuity of the molecular content found for PNe 
with very different evolutionary status and morphologies. 
In particular, N-bearing molecules are present at all the 
sampled evolutionary stages. This finding is in contrast 
with the recent suggestion made by Hasegawa et al. (2000) 
that most molecules from the AGB remnant (except H2 
and CO) are completely destroyed and then re-formed at 
the PN stage. Our observations suggest on the contrary 
that most molecular species survive during the transition 
from the AGB to the PN phase. 

An exception to this continuity in the molecular con- 
tent is provided by the sudden increase in the HCO"*" abun- 
dance. This ion is expected to form in warm PDR-like re- 
gions such as the interface between ionized and molecular 
gas. Even if currently available models fail in reproduc- 
ing the observed abundance, they clearly show that the 
abundance is sensitive to the total density (see e.g. Table 
3 in Hasegawa & Kwok 2001). Additional observational 
probes of these PDR-like regions would be provided by 
measuring the 00+ and N2H+ abundances, but the cor- 
responding observations would require higher sensitivity 
than that of the observations presented here. 



E. Josselin and R. Bachiller: The chemistry of compact planetary nebulae 



7 



4.2. Chemistry of N-bearing molecules 

The chemistry of N-bearing molecules has been recently 
studied by Ali et al. (2001) with a detailed time-dependent 
chemical moclc4 suited to the conditions of evolved PNe, 
and by Hasegawa et al. (2000) and Hasegawa & Kwok 
(2001) for the particular case of NGC 7027. The CN/HCN 
ratio calculated by Hasegawa et al. (2000) is in the range 
1.4 to 6.6, depending on the volume hydrogen density and 
the kinetic temperature, in good agreement with our av- 
erage value ('--^2.6). The abundances of HCN and HNC 
calculated by Ali et al. (2001) are in general agreement 
with the observations reported here. However their abun- 
dance of CN is too high by a factor of ^3. 

The chemistries of the three observed molecules (HCN, 
HNC and CN) are intimately interconnected. CN is mainly 
produced through the pliotodissociation of HCN or HNC, 
both having similar unshielded rates (~ 1.4 10~^ s~^, Le 
Teuff et al. 2000). It is thus expected that CN is mainly 
located in a kind of shell surrounding the HCN/HNC re- 
gions, where the ambient UV radiation coming from the 
central star can penetrate to photodissociate HCN, HNC, 
and other related compoimds. The CN and HCN emitting 
regions are thus not expected to be coincident, and the 
CN/HCN abundance ratio can strongly vary with the po- 
sition in the nebula. The CN/HCN column density ratios 
provided above must be considered as average values. 

One could ask about the potential importance of two 
body reactions in the chemistry of nitrogen. Reactions 
such as that of atomic carbon with NO and N2 or atomic 
nitrogen with CH or CS, should have negligible contri- 
butions, both because such reactions have comparatively 
lower rates and the reactants must have very low abun- 
dances. For example, atomic carbon is only present in a 
very thin transition region between molecular gas (CO) 
and ionized gas (C+), as carbon ionization potential and 
CO dissociation energy are similar (~ 11 eV). 

One exception may be the reaction 

O + HCN ^ CN + OH (2) 

as atomic oxygen may be abundant in photodissociation 
regions. Unfortunately, the abundance of atomic oxygen 
is not well known in these regions, and it will certainly 
strongly depend on the carbon-rich or oxygen-rich char- 
acter of the envelope. If atomic oxygen turned out to be 
significantly abundant, this reaction could accentuate the 
production of CN at the expense of HCN. Nevertheless, 
the presence of atomic oxygen can also contribute to the 
destruction of CN: 

+ CN ^ CO + N or NO + C (3) 

The relative importance of these processes should be ex- 
amined in the peculiar case of PN IC 5117 (see Sect. 3.2). 

In principle, HCN and HNC can convert to each other 
through hydrogen exchange reaction 

H + HNC ^ HCN + H (4) 



However, the forward reaction must be predominant since 
it has an activation energy four times smaller than the 
reverse reaction (Talbi et al. 1996). On the other hand, the 
HNC/HCN ratio is observed to decrease at increasingly 
high densities in interstellar clouds (Turner et al. 1997). 
The relatively low HNC /HCN ratio observed in PNe could 
then be due to the combination of moderate temperatures 
and densities prevailing in the molecular envelopes. 

HCN and HNC can indeed react with several molec- 
ular ions such as HCO+, Hg , H3O"'" or C2H+ to form 
HCNH+, which after dissociative recombination will lead 
to CN, HCN or HNC in a ratio 2/1/1 (Le Teuff et al. 
2000). As the abundance of molecular ions is expected to 
increase with the state of evolution (as the nebula is be- 
ing progressively ionized), such reactions are expected to 
additionally enhance the CN/HCN ratio as the evolution 
proceeds. 

4.3. Molecular ions: HCO+ and N2H+ 

As already discussed by Bachiller et al. (1997a). a remark- 
able result concerning the chemical content of PNe is the 
rather high column densities of HCO+. The "standard" 
path of HCO"*" formation which may dominate in dense 
interstellar clouds invokes the reaction of Hg with CO, 
and thus requires a high formation rate of Hg" through 
cosmic-ray ionization of H2 (Turner 1995). However the 
survival of Hg ion is expected to be very short, because of 
its very effective dissociative recombination (Bachiller et 
al. 1997a), so this process may not be efficient enough to 
explain the observed high abundances. It has been pointed 
out by Ali et al. (2001) that X-ray emission from the cen- 
tral hot star of PNe may also contribute to a high rate of 
ionization. In fact, the relatively high X-ray fiuxes have 
been observed in several PNe (e.g.: Apparao & Tarafdar 
1989) are believed to be generated from the central star 
(not from the nebula). The effect of such X-rays on the 
envelope can be particularly important in the case of com- 
pact PNe as these studied here. 

An additional way to produce HCO+ is through the 
reaction 

CO+ + H2 ^ HCO+ + H (5) 

C0+ could be formed through charge transfer between H+ 
and CO, via the reaction 

C+ + OH C0+ + H (6) 

or by the direct ionization of CO by the means, for in- 
stance, of the propagation of an ionizing front into the 
molecular gas of high CO abundance (Latter et al. 1993). 
The latter path may be dominant at least in some PNe. 
The initiation of such ionizing front is one of the charac- 
teristics of PPNe (Kwok 1993) and HCO+ seems indeed to 
form essentially during this stage (Bachiller et al. 1997a). 
Besides the formation of HCO+, the other main destruc- 
tion process of C0+ is dissociative recombination. A col- 
umn density A'^(CO^) ^ 10^^ cm~^ can be expected and 
is consistent with the detection in NGC 7027 (Latter et 



8 



E. Josselin and R. Bachiller: The chemistry of compact planetary nebulae 



al. 1993). A systematic search for this ion in PNe where 
HCO+ emission is intense would confirm this scenario. 

HCO^ is mainly destroyed by dissociative recombina- 
tion with electrons. These chemical processes are included 
in the detailed model of Ali et al. (2001) that leads to an 
abimdancc ratio HCO+/HCN^0.1, in reasonable agree- 
ment with the less evolved PNe of our sample, but sig- 
nificantly lower (by a factor of 5 to 10) that the values 
found in PNe at intermediate stage of evolution (sec also 
Bachiller et al. 1997a). The observed trend in the evo- 
lution of HCO+ abundance seems to imply an increase 
in electron density during the late stages of PN evolu- 
tion. Hasegawa & Kwok (2001) found that in NGC 7027 
about 90% of carbon, which is a major source of elec- 
trons, is in the form of C~^. Carbon with higher ioniza- 
tion degree (0^+ and 0"^+) is commonly observed toward 
more evolved PNe. An increase in the destruction pro- 
cess of HCO"'" during the last evolutionary stages is thus 
expected and consistent with our observations. The abun- 
dance of HCO+ estimated by Hasegawa et al. (2000) for 
NGC 7027 is in general agreement with our observations, 
but it assumes that it is formed in gas at a temperature 
of ~800K. Observations of spectral lines from higher J 
levels (at higher frequencies) would be useful to estimate 
the gas physical conditions. 

In molecular clouds, N2H+ is often thought to be asso- 
ciated with HCO+, since it is essentially formed through 

+ N2 ^ N2H+ + H (7) 

But contrary to HCO"'", no important alternative forma- 
tion is expected. The abundance of N2H+ was however 
not reported in the models of Ali et al. (2001). The lack of 
other significant paths of formation, together with its high 
molecular dipole moment (~ 3.4 Debye) explains why we 
only detected this ion in NGC 7027. Assuming similar line 
intensity ratios in the other PNe, one needs to reach a rms 
^ 0.5 mK to allow its detection, far below the limit we ob- 
tained. A deeper search for this ion would be interesting 
to examine any correlation with the temperature of the 
central star and in particular with its X-ray flux (when 
Z 10^ K). 

4.4. dumpiness and PDR-like regions 

Beside the evolutionary trends described above, a rather 
large dispersion of abundance ratios is observed for dif- 
ferent values of Mm/Mi. Such scatter may be, at least 
partly, attributed to the the number and density of the 
clumps, which are expected to vary from one PN to an- 
other. Indeed these variations should act on the observed 
abundances by two ways. First, the degree of dumpiness 
(filling factor) can be different for different molecules in 
the highly fragmented envelopes of the more evolved ob- 
jects such as BV 5-1 (Josselin & Bachiller 2001). Second, 
the survival of the molecules is expected to strongly de- 
pend on the number and density of the clumps which de- 
termine the shielding of dissociating UV radiation from 



the central star. Indeed, high resolution maps of the molec- 
ular envelopes of evolved PNe show that the molecular 
material is always confined within the ionized nebula (see 
e.g. Speck et al. 2002 for the Helix nebula and Josselin & 
Bachiller 2001 for BV 5-1), in contradiction with the sim- 
ple image of an initially neutral homogeneous nebula be- 
ing gradually ionized ("Stromgren sphere"). As confirmed 
by direct observations of some individual dense clumps in 
the Helix (Huggins et al. 1992, 2002), the molecules are in 
fact concentrated within dense blobs which reach a state 
of approximate dynamical equilibrium with the surround- 
ing ionized material, and whose evolution is determined 
by their density, i.e. by their capability to be efficiently 
shielded from the ambient UV radiation. Any consistent 
chemical model of the neutral matter in PNe should thus 
take into account this clumpy structure which seems to be 
an essential ingredient of PNe. 

5. Conclusions 

The presence of molecular gas in compact PNe provides a 
unique opportunity to examine a rapidly evolving chem- 
istry in a peculiar medium which is subjected to strong 
radiation fields. The observations reported in this paper 
show that some species such as HCO+ and CN are partic- 
ularly abundant compared to envelopes around AGB stars 
or even to interstellar clouds. The chemical composition 
of the molecular envelopes in these compact PNe appears 
somewhat intermediate between the composition of proto- 
PNe (such as CRL 2688 or CRL 618) and well evolved PNe 
(such as the Ring, M4-9, or the Helix). 

The molecular abundances are well explained, in gen- 
eral terms, by time-dependent gas-phase chemical models 
in which the ionization rate is enhanced by several or- 
ders of magnitude with respect to the average interstellar 
value (Ali et al. 2001, Hasegawa et al. 2000, Hasegawa 
& Kwok 2001). Thus, our observations confirm that the 
chemistry in the neutral shells of PNe is essentially gov- 
erned by the high energy radiation from the hot central 
stars. Shocks due to fast winds which develop during the 
formation of the PN can also be important, but mainly 
during the proto-PN stage and in the regions where the 
jets impact on the envelope. The complexity of the chemi- 
cal processes is increased by numerous factors linked to the 
properties of the central star and the geometry and degree 
of dumpiness of the envelope. The survival of molecules 
can only occur in well protected dense clumps, and a rich 
photo-chemistry very likely takes place at the surfaces of 
these PDR-like regions. 

Although the main processes governing the chemistry 
in the neutral envelopes of PNe seem well identified, some 
significant discrepancies with the models remain to be 
understood. For instance the calcidated CN/HCN ratio 
is a factor of ^^5 higher than the observed value, and 
the observed HCO"*" abundances tend to be higher than 
those predicted by the models. In addition, as discussed by 
Hasegawa et al. (2000), some of the molecules observed in 
evolved PNe (such as the Ring or the Helix) could not be 
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direct remnants from the AGB envelope, but newly formed 
species, this is the case of HCO+. However, the observa- 
tions suggest that there is a continuity in the chemical 
composition from the proto-PN phase to the full photodis- 
sociation of the envelope, and that molecules like CO and 
HCN arc progressively photo-dissociated. 

Additional observations are now required to push up 
our understanding of this chemistry. The search for species 
such as CO"'" woiild be useful to confirm the chemical net- 
works invoked to explain the observed abundances. 
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